Excitation energy dependence of symmetry energy of finite nuclei 
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Abstract 

A finite range density and momentum dependent effective interaction is used to calculate the 
density and temperature dependence of the symmetry energy coefficient C sym (p,T) of infinite nu- 
clear matter. This symmetry energy is then used in the local density approximation to evaluate the 
excitation energy dependence of the symmetry energy coefficient of finite nuclei in a microcanonical 
formulation that accounts for thermal and expansion effects. The results are in good harmony with 
the recently reported experimental data from energetic nucleus-nucleus collisions. 
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The symmetry energy E sym (p,T) represents with a very good accuracy the energy cost 
per nucleon to convert all the protons to neutrons in symmetric nuclear matter at the 
density p and temperature T. Study of its density and energy dependence is of utmost 
contemporary importance. It is essential not only for understanding many aspects of exotic 
nuclear physics induced by collisions of radioactive nuclei, but also a number of important 
issues in the astrophysical scenario like supernovae explosions [l[ , explosive nucleosynthesis, 
cooling of protoneutron stars and abundances of relatively heavier elements. Even on a 
more mundane level, the neutron skin thickness of heavier nuclei is intimately correlated to 
the density derivative of the symmetry energy 0, 0] as it reflects the pressure difference on 
the neutrons and protons. 

In addition to a kinetic contribution, the symmetry energy has a contribution arising from 
the difference between the neutron-proton (n-p) interaction and that between like pairs (n-n, 
p-p). Given an interaction, it is straightforward to calculate the symmetry energy at different 
densities and temperatures for infinite matter. There have been several attempts in this 
direction. Calculations of the nuclear equation of state (EoS) in the microscopic framework 
using both bare [3, lEf and effective interactions 0, 0] have been done. The outcome of these 
calculations for the symmetry energy is similar (~ 30 — 35 MeV) at saturation density, but is 
considerably different at subnormal as well as at supranormal densities where the available 
data from experiment to confront with theory are more scarce. 

Laboratory information on the density dependence of the symmetry energy can be ob- 
tained from energetic nucleus-nucleus collision experiments. At densities above the normal, 
it can be inferred from the comparison of theoretical predictions with experimental data on 
the differential flow of neutrons and protons, from the -k~/tt + , K°/K + ratios, etc 0. At 
subnormal densities, disassembly of a hot expanded nucleus offers the best tool to study 
the characteristics of the symmetry energy p, T(|. Experimental data related to isotopic 



distributions, isospin diffusion, and isoscaling try to constrain the density dependence in the 
subnormal region, but there is still considerable uncertainty. 

A nucleus expands with excitation with increasing temperature in general. This implies 
an excitation energy dependence of the symmetry energy because of the density change. 



Experimentally, this information is generally extracted [lOj from the fit of the experimental 
isotopic distributions at different excitation energies to those obtained from a model for mul- 
tifragmentation like the statistical multifragmentation model (SMM) [ll[ or from isoscaling 
12| . Currently, calculations for the energy dependence of the symmetry energy are available 
for infinite matter, but no microscopic calculation has yet been performed for the energy 
dependence of finite nuclei. The main purpose of this communication is to report such a 
calculation. 

For an expanding system pursuing the equilibrium configuration, the surface diffuseness 
is likely to play an important role [Lj], thus a zero-range interaction like the Skyrme force 
widely used to explore the nuclear ground-state properties may not be the most adequate for 
generating such a density profile. It is further noted that a constrained expanded system in a 
Thomas-Fermi approach may lead to numerical instabilities [14] and the gradient (surface) 
lerms in the energy density functional were replaced with a suitable Yukawa interaction 
15|. We have therefore chosen the modified Seyler-Blanchard (SBM) effective interaction 
16j for our microscopic calculation in the finite temperature Thomas-Fermi formulation. 
This interaction is of finite range and momentum and density dependent. The interaction 
reproduces quite satisfactorily the ground state bulk properties of nuclei over the whole 
periodic table for A > 16. The EoS calculated 17] with this interaction agrees very favorably 
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with those obtained microscopically with a realistic interaction in a variational approach 
18l . Il9l |. The SBM interaction is given by 



v(r 1 ,r 2 ,p,p) = -Q )U 



l-^-rf 2 {p(ri)+p(r 2 )}" 



exp(— r/a) 
(r/a) 



(1) 



Here r = |ri — r2 1 and p = |pi — P2I are the relative separation of the interacting nucleons 
in coordinate and momentum space, p(ri) and p(r 2 ) are the densities at the sites of the two 
interacting nucleons, and C\ and C u are the strengths for like pair and unlike pair nucleon- 
nucleon interaction. The density exponent n controls the stiffness of the nuclear EoS. The 
values of the parameters Cy, C u , b, d, a and n are given in Refs. [l6|, 2(j. The energy per 
nucleon e(p, T) calculated with this interaction for infinite nuclear matter is given by 

e( P ,T) = -^Pr [TJs/2for)M/2for) {l ~ ^X} + \v?\ , (2) 

P T L J 

where r refers to the isospin index (n, p). Here J q (r)) are the Fermi integrals, m^. the effective 
A;— mass of the nucleon and i] T the fugacity given by 



Vr = (fir ~ K° - V?)/T, 



(3) 



with p T as the nucleon chemical potential. In Eqs. (2) and (3), the VVs are the different 



components of the single-particle potential whose expressions can be found in [16 



The symmetry energy per nucleon of asymmetric nuclear matter with asymmetry X 

(Pn -Pp)/P IS 



It can be written as 



e sym (p, T, X) = e(p, T, X) - e(p, T, X = 0). 



e sym (p, T, X) = C sym (p, T)X 2 + 0(X 4 



(4) 



(5) 



The terms beyond X 2 are negligible over a considerable range of X (as involved in finite 
nuclei). The symmetry energy coefficient C sym is obtained from 



1 d 2 

C sym (p,T) = --g-^e sym (p, T, X)\ x =o- 



(6) 



In the top panel of Fig. 1, the density dependence of C sym at T = is displayed in 
the density region 0.1 < p/po < 1.0 for three variants of the SBM interaction with density 
exponents n = 1/6, 2/3, and 4/3, in increasing order of the stiffness of the nuclear EoS. 
The values of the nuclear incompressibility with these three interactions are = 238, 300, 
and 380 MeV. The symmetry energy coefficient calculated with them can be very nicely 
represented by C sym (p) ~ C sym (p ) (p/p ) 7 with C sym (p ) = 34.0 MeV and 7= 0.65, 0.68, 
and 0.70, respectively. The value of the exponent 7 appears not very sensitive to the nuclear 
EoS. The agreement of the functional form of the symmetry energy coefficient and the value 
of 7 with those obtained recently [l(| from experimental data (7 ~ 0.69) is excellent. The 
nuclear incompressibility with n — 1/6 compares very well with the presently accepted 
value of ~ 230 MeV; all the subsequent calculations are therefore reported for the SBM 
interaction with n — 1/6. In the bottom panel of Fig. 1, the temperature dependence of the 
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symmetry energy coefficient at different densities (p/po = 0.3, 0.6, and 1.0) is shown. At a 
fixed density, dependence on temperature is not much evident. 

In calculating the excitation energy or density dependence of the symmetry energy co- 
efficient of a finite nucleus with N neutrons and Z protons (A = N + Z) with excitation 
energy E*, we remind ourselves that the hot nucleus prepared in the laboratory in energetic 
nuclear collisions is an isolated system with a fixed total excitation E* and thus should be 
described by micro canonical thermodynamics [2lj]. Left to itself, the system expands due 
to unbalanced thermal pressure in search of maximal entropy where the total pressure van- 
ishes and the system is in equilibrium in a bloated mononuclear configuration with the same 
excitation energy. The expansion is simulated through a self-similar scaling approximation 
for the density: 

p x (r) = A 3 p(Ar), (7) 

where A is the scaling parameter (0 < A < 1) and p(r) is the base density profile. The base 
density, employing the SBM interaction, is generated in the self-consistent Thomas-Fermi 
framework. The subtraction scheme [22, 23[ is used to render the density profile independent 



of the box size with an effective temperature T chosen so as to give the maximum entropy 
for the given excitation E*. The excitation energy is calculated as 

E* =E(\ eq ,T)-E(\ = l,T = 0), (8) 

where X eq is the scaling parameter for the equilibrium density profile corresponding to this 
excitation. 

The SBM interaction, being momentum dependent, renormalizes the bare nucleon mass 
m to an effective k— mass. A frequency dependent mass factor m^/m is further phenomeno- 



logically incorporated [241 . 125] in the calculation. It is very relevant in the present context; 
the o>-mass m^/m is generally larger than unity, it has the effect of bringing down the excited 
states from higher to lower energy near the Fermi surface, thus increasing the many-body 
density of states at low excitations that allow comparatively more accommodation of entropy 
at a given excitation energy. Details on the generation of the equilibrium density profile, 
effective temperature, frequency dependent mass, etc., as employed in this calculation, are 
given in Refs. 0, 2?] . 



Once the equilibrium density p(r) of a nucleus at excitation E* is known, the symmetry 
energy is calculated in the local density approximation as 

CUE') (^) 2 = i / P(r) C M ,T) (e&Z^y*. (9) 

Here C l sym (p(r),T) is the symmetry energy coefficient at temperature T of infinite nuclear 
matter at a value of the local density p(r). The local isospin density is given by p n {f) ~Pp{ r )- 
It may be mentioned that both the volume and the surface terms in the liquid drop type mass 
formula are asymmetry dependent |28|. The symmetry energy coefficient C sym (E*) defined 
through Eq. (9) may therefore be taken as an effective parameter incorporating both the 
volume and surface contributions from asymmetry and may be written as 

/~i _ fivol _ (~<surf I aI/3 /-irv\ 
^sym — ^ S ym ° sym I • \ LU ) 

In a microcanonical formulation, it has been found that the equilibrium density at a 



given excitation depends on the mass and asymmetry of the nucleus concerned 26]. In 
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investigating the excitation energy dependence of the symmetry energy coefficient, it would 
then be worthwhile to investigate its system dependence. We have therefore chosen three 
systems, two isobars of A = 150, namely, Cs and Sm, and a lighter system 40 S. In the 
bottom panel of Fig. 2, the coefficient C sym as calculated using Eq. (9) is displayed as a 
function of E* /A for the isobars of A = 150. It is found that at a fixed excitation, including 
the ground state, C sym is somewhat sensitive to the asymmetry of the nucleus; it increases 
with increasing proton fraction of the system. This is at variance with the expectation from 
the liquid-drop formula where the effective symmetry energy coefficient given by Eq. (10) 
is independent of charge for a given mass. This may be understood from the fact that the 
parameters of the liquid-drop formula are based on a global fit to the binding energies of 
nuclei over the entire periodic table around the stability line excluding the very light ones. 
Here, the Coulomb energy (= a c Z 2 /A 1 ^ 3 ) coefficient a c is taken as a constant for the whole 
mass range. In our calculations for isobars, it is seen that with increasing charge as the 
proton distribution is pushed outward, a c decreases and C sym increases. This effectively 
explains the isobaric variation of binding ener gies . 

Some representative experimental results [lOj for C sym obtained from the analysis of 
isoscaling data for lighter fragments are shown in the bottom panel of Fig. 2 as inverted 
open triangles and solid circles. At relatively lower excitations (E*/A ~ 2-3 MeV), ex- 



perimental data from isoscaling for heavier fragments [29] are also shown as solid squares. 
Our calculated energy dependence of the symmetry energy coefficient agrees favorably with 
these experimental findings, the calculated values are somewhat lower. The variation of C sym 
with excitation stems basically from the changing equilibrium density with excitation en- 
ergy, which is shown in the middle panel of the figure for the two systems. Since the density 
has a profile, the choice of a single value of the density leaves room for ambiguity; we have 
taken p c j p Cj o as the measure of the density where p c is the central equilibrium density at the 
relevant excitation and p Cj o is the ground state central density. The theoretically calculated 
results are in nice agreement with the experimental data from Ref. [3(| (solid circles of the 
middle panel) derived from the analysis of caloric curve measurements. The data obtained 
from Coulomb barrier systematics 3l|, |32) are shown with open squares. For completeness, 
in the top panel, the caloric curves along with the experimental data compiled by Cibor et al 



33J are also displayed. It is seen that the plateau of the caloric curve shows little sensitivity 



to the asymmetry of the nucleus; this is consistent with the recent calculations of Hoel et al 
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The mass dependence of the effective symmetry energy C sym (E*/A) is displayed in the 
bottom panel of Fig. 3. The calculations are done for 150 Sm and 40 S, both having nearly 
the same asymmetry. The reduction in C sym for 40 S can be understood from the role played 
by the surface asymmetry as given by Eq. (10). For completeness, the equilibrium central 
densities and temperatures as a function of E*/A are also shown in the middle and top panel 
of the figure, respectively. 

To conclude, calculations on the density and excitation energy dependence of the sym- 
metry energy of finite nuclei have been reported in this communication in a microscopic 
formulation within the microcanonical framework. It has been stressed in a recent calcula- 
tion |l3j that the surface diffuseness of the expanded mononuclear system plays a key role 
in making the system softer towards instability, limiting the maximum excitation energy a 
mononucleus can hold to ~ 5 MeV/ A with free variation of the surface diffuseness. Our 
model calculation does not leave any room for free variation of the surface diffuseness. It 
is determined in two stages: the increased diffuseness of the base density profile of the hot 
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nucleus over that of the ground state and then its subsequent stretching from the self-similar 
expansion. The surface diffuseness so obtained is found to be somewhat less than that re- 
ported in 13, 34j. Exploring the weakening of the symmetry energy with excitation using 
free variation of the surface diffuseness would be interesting to look into. In the present 
calculation, the density dependence of C sym (p) of infinite nuclear matter is found out to be 
~_lp/Po) 7 with 7 = 0.65, very close to the recently extracted experimental value of 7 ~ 0.69 
[lOj. At constant density, C sym (E*) of infinite nuclear matter is practically constant. For 
finite nuclei, however, density changes with excitation; their excitation energy dependence 
can be well represented by C sym (E*) ~ C sym (E* = 0)(1 - aE*) with a ~ 0.06. These are 
in good consonance with the experimental data obtained from nuclear multifragmentation. 
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Figure Captions 



Fig. 1 Symmetry energy coefficient of infinite nuclear matter as a function of density for 
different variants of the SBM interaction (see text) at T = 0.0 MeV (top panel). In 
the bottom panel the temperature dependence of C sym at several fixed densities is 
shown for n = 1/6. 

Fig. 2 The equilibrium temperature (top panel), equilibrium central density (middle panel) 
and the symmetry energy coefficient (bottom panel) as a function of excitation energ 
for the A = 150 isobars (Cs and Sm). The experimental data for T are from Ref. [3c 



those for p/p are from Refs. [30J (circles) and [3l|, [32| (squares) , and those for C, 



are from Refs. [lOj (inverted triangles and circles) and [29] (squares). 
Fig. 3 Same as in Fig. 2 for the systems 40 S and 150 Sm to show the mass dependence. 
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